m 



$ w NNY VA t E, C AL1 F t A 


N73-31638 

Unclas 

15396 


/£LflSC~D350709) ANALYSIS OF LOft ENERGY 
ELECTRONS Final Report {Lockheed 
Hissiles and Space Co* J ^ 15 p HC $3 O Q0 

CSCL 20 fl 

! G3/24 


Enclosure (l) to LMSC/D 350709 
1 June 1973 


FINAL EE PORT 

ANALYSIS OF LOW-ENERGY ELECTRONS 


Contract NASw 2253 


Prepared “by: 


Dr. Richard D_. Sharp 
Space Sciences Laboratory 
Lockheed Palo Alto Research Laboratory 
Lockheed Missiles & Space Company, Inc. 
3251 Hanover Street 
Palo Alto, California 9^304 


x 



TABLE OF CONTENTS 


Section Title Page 

1 Lockheed Experiment Description. • 1 

2 Computer Programs 2 

2.1 Survey Program 2 

2.2 Plasma Properties Program 4 

3 Analytic Approximations. 11 

4 Results 11 


i 



F INAL REPORT 

ANALYSIS OF LOW-ENERGY ELECTRONS 
Contract NASw 2253 


This was a cooperative research program undertaken by the Lockheed 
Palo Alto Research Laboratory (LPARL) and the Massachusetts Institute of 
.Technology- (MIT ) to -analyze- simultaneous observations of low-energy elec- 
trons in the plasma sheet and in the auroral zone. Data from the MIT 
Plasma Experiment on the OGO-3 satellite and from the Lockheed experiment 
on the 0V1-18 satellite ( 1969 - 25 B) were processed and intercompared. The 
principals involved were Professor Vytenis M. Vasyliunas of MIT and Dr. 
Richard D. Sharp of Lockheed. This report describes the effort expended 
and the results obtained under the Lockheed portion of the program. 


1. Lockheed Experiment Description 

The Air Force satellite 0V1-18 was launched on l8 March 19&9 i^to a 
99° inclination orbit with apogee at 590 km and perigee at U 69 km. An 
experiment by the Lockheed Palo Alto Research Laboratory included an assort- 
ment of low-energy particle detectors designed to measure the intensity, 
angular and energy distributions of the auroral electrons and protons in 
the energy range below about 50 keV. The detectors of primary interest 
for the comparison with the MIT plasma experiment are a set of 13 magnetic 
electron spectrometers. They utilize channel-electron -multipliers as the 
sensing elements and cover the energy range from 0.8 to 37 keV. They are 
similar in design to the Lockheed detectors flown on the ATS -5 satellite 
and have been described in the literature . A matched set of four of 
these detectors was positioned at each of three angles spanning a range 
of 90°. The characteristics of the individual detectors are shown in 


"A Low-Energy Channel-Multiplier Spectrometer for ATS-E," R. D. Reed 
E. G. Shelley, J. C. Bakke, T. C. Sanders and J. D. McDaniel, IEEE 
Trans. Nucl. Sci ., NS-l6 , 359? Feb. 1969 ; "Channel -Multiplier Instru- 
mentation for the Measurement of Low-Energy Auroral Particles," M. F. 
Shea, G. B. Shook, J. B. Reagan, L. F. Smith, and T. C. Sanders, IEEE 
Trans. Nucl. Sci ., NS-lU , 96, 1967. 
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Table I. The other instruments comprising the satellite payload included 
a set of three scintillator -photomultiplier threshold detectors which meas- 
ured the lover energy electron fluxes (thresholds at 0.2 and 1.2 keV) at 
three different angles, and various foil threshold detectors measuring the 
more energetic electrons as well as an assortment of proton detectors. Due 
to a malfunction in the stabilization system, the satellite is slowly tumbling. 
A three-axis magnetometer provides aspect information, however, and for the 
purposes of this intercomparison, no serious degradation of the data results 
from the tumbling. 

2, Computer Programs 

A substantial amount of effort under this contract was devoted to the 
development of two computer programs which were required for the reduction 
of the 0V1-18 data. These are described below. 

2.1 Survey Program 

Modifications and revisions were made to a previously- existing computer 
program for unpacking the Air Force- supplied digital tapes and producing 
survey plots showing the outputs of the various detectors on the satellite 
as a function of time on a scale of approximately 100 seconds per cm. These 
plots are needed to evaluate the quality of the data and the magnitude of 
any background corrections which are applied in later stages of the analysis. 
The program also produces as output a computer-compatible digital tape con- 
taining the satellite data in a suitable form for further analysis on the 
Lockheed computers. 

As has been indicated, the 0V1-18 satellite was tumbling and this 
introduced occasional backgrounds to the detectors due to scattered ultra- 
violet light and to the ambient cold plasma when the detectors were point- 
ing in the direction of the velocity vector. In order to evaluate these 
effects, the outputs of all 48 low-energy particle detectors contained in 
the payload, were included in the survey plot outputs. These included 
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Table I. CHARACTERISTICS OF PERMANENT MAGNET- ELECTRON SPECTROMETERS ON OV1-18 


Detector 

Energy- 

Range 

Geometric 

Factor 

Nominal 

Pointing 

Angle 

(With 

. Respect to 
Zenith) 

\ 

Angular 
Acceptance 
Full Angle 
in Degrees 

Sensitivity 

Threshold 

Particles/ 

cm^-sec 

■ Group 1 
CME-1A 

0.8-1. 5 

0.79x10"^ 

0° 

10x10 

4x10;? 

CME-1B 

1.75-3.3 

0.8lxl0“' 

■ 0° 

10x10 

3x10;? 

CME-1C 

3.75-7.0 

0 . 96 x 10 "? 

0° 

r\ 

10x10 

3x10' 

CME-1D 

8 . 3 - 16.3 

0 . 26 x 10 7 

0° 

10x10 

1x10? 

CME-1E . 

17.3-37.0 

0.40x10 

0° 

0 

10x10 

5x10 

CME-1F 

Penetrating 

— 

0 

— 

--** 

Group 2 
CME-2A 

Background 
0.8-1. 5 

0.67x10"^ 

55° 

10x10 

3x10;? 

CME-3B 

1 . 8 - 3. 3 

0.84xl0~j? 

55 o 

10x10 

3xio^ 

CME-3C 

3.75-7.0 

0 . 98 x 10 ? 

55° 

10x10 

2x10? 

1x10^ 

CME-2D 

8 . 2 - 16.5 

0.22x10 

55 0 

10x10 

CZU-2 

Ultraviolet 

— 

55 

— 

— 

Group 4 
CME-4A 

Background 

0 . 8 - 1.5 

0 . 69 x 10 '^ 

90 ° 

10x10 

4x10;? 

CME-4B 

1 . 8 - 3.25 

0.94x10 ' 

9° 

10x10 

? xl0 5 

CME-4C 

3. 7-7.0 

0 . 85 x 10 ? 

90° 

10x10 

4x10? 

CME-4D 

8.4-16.3 

0 . 25 x 10 4 

90° 

10x10 

7x10 




ultraviolet and penetrating radiation background detectors. A typical 
section of this output (for one of the coordinated cases), including data 
from four of the electron spectrometers used in this analysis, is shown 
in Figure 1. The voltage outputs shown are approximately logarithmically 
related to count rate over a four-decade dynamic range. Ephemeris inform- 
ation, magnetic local. time (AT = auroral time), and L values are computed in 
a separate program. 

2.2 Plasma Properties Program 

A computer program which takes the unpacked digital tape produced as 
output by the survey program and computes the various plasma properties at 
one-second intervals as a function of Universal Time (UT) and pitch angle 
was developed specifically for this project in a form agreed upon in the 
course of discussions with Vasyliunas. This is a versatile program which 
has provisions for easy modification of spectral shape assumptions, back- 
ground levels, detection efficiencies, etc. It produces a series of ten 
■plots ■‘showing, -as a function of time: 

• The gross counting rates of the thirteen relevant detectors 
(3 plots); 

• Pitch-angle distributions at one-second intervals (l plot); 

® Electron spectra at each of three pitch angles at one- 

second intervals (3 plots); and 

• The electron plasma properties: number flux, number dens- 

ity, flux -weighted average energy and density-weighted 
average energy, at each of three pitch angles at one- 
second intervals (3 plots). 

In addition, it sorts the events into 15 categories depending on the 
number of channels responding above their in-flight calibration levels in 
order to allow for the utilization of different approximations for the 
different classes of events. Some typical outputs of this program are 
shorn in Figures 2 through 5 - 
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I* 2.2 ,2.9 ^.2 6.9 IB. 5 37.3 244. 4 178.7 34.8 i4.2 8.0 

AT(HR) 23.16 23.11 23.03 22.91 22.69 22.19 20.30 14.26 12.61 12.15 11.94 


FIGURE 1 
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Figure 2 shows the count rates of the group 2 electron detectors for 
a portion of the pass illustrated in Figure 1. , They are semi-log plots 
with the ordinate beginning at one count/sec. Each division represents 
one decade and the curves are offset from each other by two decades. The 
abscissa is UT in seconds. 

Figure 3 shows the pitch-angle distributions of the integral number 
fluxes in the range (0.8 ^ E ^ 17 keV). The asterisks represent the group 
1 detectors, the circles the group 2, and the crosses the group 4 (See 
Table l) . Each square represents a one- second measurement with ten of 
them stacked vertically in a region of abscissa corresponding to the ten 
seconds of UT over which the measurements were acquired (time progresses 
from top to bottom). The UT scale matches that of Figures 2, 4 and 
The abscissa of each individual square is pitch angle over the range 0° 
(particles coming down the field line) to l 80 ° (particles going up). 

Since the vehicle is tumbling these values change slowly with time. The 
ordinate is flux on a logarithmic scale from 10^ to 10 10 (cm 2 -sec-sterad)“ 1 

Figure 4 shows normalized electron spectrums from the four group 2 
detectors with a stacking arrangement similar to Figure 3. Energy is 
plotted approximately logarithmically along the abscissa of the individual 
squares (See Table I for the actual energy values). The ordinate is a 

two-decade logarithmic scale normalized to the peak value in units of 

/ 2 .-1 
(cm -sec-sterad-keV) 

Figure 5 shows the computed plasma properties from the group 2 detect- 
ors (0.8 ^ E ^ 17 keV) as a function of UT. ■ From top to bottom are: l) 

integral number flux on a four-decade log scale over the range 10^ to lO 1 ^ 

f 2 s —l V 

(cm -sec-sterad) ; 2) number density on a four-decade log scale over the 
range 10 to 10 (cm ); 3) flux -weighted average energy on a linear 
scale from 0 to 20 keV; 4) density weighted average energy on the same 
scale as 3); and 5) pitch angle on a linear scade from 0° (particles coming 
down) to l8o° (particles going up). 
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FIGURE 2 
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FIGURE 3 
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0V1-18 GRP 2 ELECTRON SPECTRA 


A 3 1 0 



FIGURE 4 
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3. Analytic Approximations 


Several different approximations were investigated for the computation 
of the plasma properties from the outputs of the four or five broadband 
spectrometers oriented at any one pitch angle. Initial calculations were 
done with an equivalent rectangular response function for the spectrometers 
and a simple histogram approximation to the electron spectrum shape. The 
upper channels were fairly broad, however, particularly in the group 1 detect- 
ors (See Table i) and since the flux was generally falling off rapidly with 
increasing energy in this range, we investigated some more realistic approx- 
imations. In one study, an assumed exponential flux shape was fitted to 
the two upper channels of the group 1 detectors by using the laboratory- 
measured response functions and the plasma properties were computed from 
this exponentially falling flux at the higher energies and the assumed his- 
togram flux shape at the lower energies. We processed considerable data 
to compare this approximation with the more simple approximation and con- 
cluded that the differences were generally in the range of 20 to 30 percent. 
We felt that this precision was sufficient for satisfying the principal 
objectives of this program and since there was a decided advantage to using 
the same approximations for all angular groups, we settled on the histogram 
approximation for the final analysis. 

4- Results 

A preliminary field line mapping technique was developed by Vasyliunas 
in order to project the OGO-3 position down to the altitude of the 0V1-18 
orbit and determine the approximate coordination time. We utilized this 
technique to select the most favorable half of the initially available cases 
(when OGO-3 was located in. the plasma sheet at local times near those at 
which 0V1-18 crossed the auroral zone) for processing and examined the 
survey plots in each case in order to determine the quality of the data, 
the approximate flux intensity, the degree of spatial and teuporal struc- 
ture and the magnitude of the background phenomena. Then we processed 
these cases through our plasma properties program. At that point (Dec. 1971) 
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we decided to postpone any further effort until the OGO data became avail- 
able. -Due to a changeover, in computers at MIT, -there was an unavoidable 
delay of about one year before these data could be processed. At that 
time (Dec. 1972) it became apparent that on more than half of the original 
cases, OGO data would not be available because of data gaps, missing tapes, 
etc. A total of 32 cases remained and we processed the 0V1-18 data for 
those of them that had not previously been analyzed and did the further 
detailed analysis of the in-flight calibration data necessary to insure 
the quantitative precision of the computed plasma properties. 

At that time (April 1973) we exchanged data with Vasyliunas and were 
able to make our first preliminary intercomparisons. The results look 
quite promising. Clearly apparent signatures of the expanding and contract- 
ing plasma sheet are evident in the OGO data. The 0V1-18 data generally 
exhibit the high degree of variability associated with substorm controlled 
phenomena. After discussion with Vasyliunas we agreed that the next step 
should be a study of ground-based magnetograms during the period of interest, 
and a further refinement of the field line tracing technique in order to 
establish more accurately the time of the coordination. Vasyliunas has had 
considerable experience in both these areas and he indicated that he would 
undertake these tasks. At this point our resources under this grant are 
expended and we cannot undertake any further extensive analysis effort. 
However, we will participate in the final writeup of the results for publi- 
cation. 


12 



